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Abstract 


Despite being an essential micronutrient, copper is also a potentially toxic heavy 
metal. Using selection experiments, we produced Chlamydomonas reinhardtii 
populations with increased tolerance of copper ions and then derived pure cell 
lines from these populations. Strains derived from the same population (both 
adapted and nonadapted) significantly differed in terms of growth parameters. 
Cultivation of the strains in a range of copper ion concentrations revealed 
differences in growth and photosynthetic performance, which could be attributed 
to microevolutionary processes occurring with each cell division. Our results 
demonstrate the effects of environmental factors on rapidly multiplying 
microorganisms. 


Keywords 
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1. Introduction 


Copper is an important microelement present in many redox-active proteins, 
including electron transport chain components, antioxidant enzymes, and oxidases, 
among others (Yruela, 2009). High redox activity leading to the generation of 
reactive oxygen species (ROS) is the primary reason for the toxicity of high copper 
concentrations in plant cells (Szivak et al., 2009). High intracellular copper 
concentrations induce changes in gene expression, antioxidant enzyme activity, and 
small-molecule antioxidant synthesis (Burda et al., 2003, 2006; Jiang et al., 2016; 
Sabatini et al., 2009); therefore, copper homeostasis must be stringently regulated 
through a complex system. Indeed, intracellular copper ion concentration is 
precisely regulated by specific transporters (ATP-binding cassette transporters and 
copper-transporting P-type ATPases) and chaperones (yeast Atx1 homologs) 
(Burkhead et al., 2009), and all components of this mechanism are regulated by 
several genes. Notably, this mechanism may be modified during the process of 
adaptation. 


Chlamydomonas reinhardtii P. A. Dang. (Volvocales: Chlorophyceae) can adapt to 
different copper ion concentrations. In a previous study, we showed that following 
the long-term cultivation of Chlamydomonas reinhardtii cells in a medium 
containing high copper concentrations, their ability to tolerate this metal increased 
significantly. However, the microevolutionary processes underlying this adaptation 
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slightly differed among the cultivated populations (Plucinski et al., 2018). Therefore, 
the next step was to investigate the response of cell strains obtained from the same 
population to excess copper ions. For this, we derived strains from single cells using 
the method of serial dilution. In the present study, we compared cell strains obtained 
from populations adapted to high copper concentrations (Cu) and from populations 
without previous exposure to toxic copper ion concentrations (N). In these strains, 
we studied photosynthetic performance and growth rate during cell division in 
media supplemented with different copper ion concentrations based on chlorophyll 
fluorescence measurements. The observed differences indicated that 
Chlamydomonas reinhardtii cells present various adaptations, and in selected 
populations, many combinations of microevolutionary changes could be identified. 


2. Material and Methods 
2.1. Growth Conditions, Strain Selection, and Monitoring 


For the experiments, we used the cell wall-less mutant (CW 15) of Chlamydomonas 
reinhardtii and aseptically cultured it in Erlenmeyer flasks (250 mL) containing the 
Sager-Granick medium (Harris, 2009) supplemented with 100 mM mannitol as an 
osmoprotectant and sodium acetate (75 mM) and citrate (1.7 mM) as organic carbon 
sources. The copper ion concentration of the medium was 0.25 uM. The cultures 
were incubated in a growth chamber with gentle shaking at 22 + 2 °C under the 16:8 
light:dark cycle (at 50 pmol m~? s~! photosynthetic active radiation from a 
fluorescence lamp). Cell growth was monitored by measuring the maximal 
chlorophyll fluorescence intensity, as previously described by Pluciriski et al. (2018). 


After 4 years of cultivation at high copper ion concentration (5.25 uM), from the 
population adapted to copper, we obtained strains Cul and Cu2 using the serial 
dilution method. Using the same method, we also obtained strains N1 and N2 from 
the nonadapted populations of Chlamydomonas reinhardtii. The introductory 
characteristics of these basic populations have been reported by Plucirski et al. 
(2018). The microcultures used for testing the tolerance of the selected strains to 
increasing copper ion concentrations (0.25-200 mM) were aseptically grown in 
standard 48-well flat-bottom plates without shaking. Chlorophyll fluorescence 
parameters [F,, - maximum chlorophyll fluorescence and F\,/F), - maximum 
quantum yield of photosystem II (PSIT)] were recorded using the PAM 101 (Walz, 
Germany) fluorometer. Fluorescence measurements were obtained from four 
biological replicates. 


2.2. Growth Curve Analysis 


To analyze the changes in FM in the subsequent days of the experiment, we used a 
logistic growth model (Figure 1). This model is the simplest approach to simulate 
population growth in environments of limited carrying capacity (e.g., Gotelli, 2008). 
We used the following equation, where Fy is expressed as a function of time [F),(t)]: 


Fy,(t) 


7 Asymptote 
~ 1+ exp [—Initial Growth (t — Inflection Point)|’ 





The above equation includes three parameters: 


+ Initial growth (day~’): Intrinsic growth rate (rate of increase in chlorophyll 
fluorescence); 

« Asymptote [F,, (a.u.)]: The carrying capacity of the environment and/or the 
ability of the alga to use resources; 

¢ Inflection point (day): The time point at which the growth rate of the population 
is the fastest, following which it slows. 


The above parameters of logistic growth were estimated for each strain following 
exposure to a range of copper ion concentrations. At each copper ion concentration, 
the parameters of the strains and experimental populations (adapted and 
nonadapted) were compared. The analysis of parameter means allowed for the 
identification of the strains or populations that significantly differed from one 
another. 
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Figure 1 Schematic drawing of the logistic population growth rate, with marked 
parameters estimated by the model and used in the present study. 


2.3. Photosynthetic Apparatus Efficiency Analysis 


To analyze the efficiency of the photosynthetic apparatus (F\/F,,), we used data from 
fourth day onward, when the chlorophyll levels in the samples allowed for the 
measurement of fluorescence. The data were analyzed using a linear model, which 
considered the experimental population [Cu - copper-adapted and N - nonadapted 
(control)], copper ion concentration, and time. In addition, the effects of the strains 
were nested within the population from which they were derived. Initially, all 
possible interactions were included in the model, but nonsignificant interactions 
were sequentially removed to increase the power of the test. The parameters fitted to 
the linear model are listed in Table 1, Table 2, and Table S2. 


2.4. Software 


Fy and F\,/F,, were analyzed using JMP 12.2. All plots were created using the 
R package ggplot2 — a system for declaratively creating graphics (Wickham, 2016). 


3. Results 
310 Fe 


FM measurements revealed differences in the three parameters of population growth 
(for raw data, see Figure 2). The most pronounced differences were observed in 
growth asymptote, which differed between the Cu and N populations at all copper 
concentrations. The Cu populations always exhibited higher asymptotic growth than 
the N populations, and the magnitude of these differences increased with increasing 
copper ion concentration (Figure 3) (see Table S1 for parameter estimates). At the 
three highest copper ion concentrations, the strains differed in terms of asymptotic 
growth, and these differences were more pronounced in the Cu populations than the 
N populations (Figure 3). 


The initial growth of the N populations was affected by all copper ion 
concentrations. However, in the Cu populations, decreased initial growth was 
detected only at the three highest concentrations (see Figure S1 for detailed graphs). 
The initial growth rate of strain N1 was significantly lower than that of the other 
strains at copper ion concentrations of 25 (1.22 day!) and 50 uM (1.06 day“). 
Meanwhile, the initial growth rate of strain N2 was lower at copper ion 
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Table 1 Analysis of the maximum quantum yield of photosystem II (Fy/F,) considering 
the effects of the strains nested within the experimental population, measurement day, and 
copper concentration, as well as the relevant interactions of these factors. 





af F Pp 

Experimental population 1 12.88 0.0004 
Strain (experimental population) 2 17.45 <0.0001 
Day 1 73.71 <0.0001 
Copper 1 436.90 <0.0001 
Experimental Population x Day 1 10.25 0.0015 
Day x Copper 1 7.40 0.0068 
Error 376 


df — degrees of freedom; F - calculated statistics; p - probability of obtaining a result by chance. 


Table 2 Separate analyses of the maximum quantum yield of photosystem II (Fy/F),) in 
the two experimental populations. 








Cu populations N populations 

af F p af F Pp 
Day 1 92.45 <0.0001 1 11.70 0.0008 
Copper 1 253.85 <0.0001 1 200.25 <0.0001 
Day x Copper 1 11.31 0.0009 1 0.70 0.4033 
Strain 1 35.20 <0.0001 1 6.83 0.0097 
Error 187 187 


In each case, two strains were exposed to a range of copper concentrations and their Fy/F,, values were 
measured on subsequent days. df - degrees of freedom; F - calculated statistics; p - probability of 
obtaining a result by chance. 


concentrations of 75 (1.15 day’) and 100 uM (1.06 day~'). At the highest copper 
ion concentration, the N populations did not grow at all, and the two Cu strains 
differed in terms of the initial growth rate (Cul: 1.36 and Cu2: 0.85 day‘; Table $1). 
Specifically, strain Cul grew 60% faster than Cu2. 


The time point at which the populations reached the inflection point (i-e., exhibited 
the fastest growth, just before slowing down) fluctuated with copper ion 
concentration. At copper ion concentration of 25 uM, the growth rate of both N 
strains continued to accelerate for a relatively long time (N1: 4.78 and N2: 4.75 days). 
However, at copper ion concentrations of 100 and 125 uM, this duration was 
shortened for the N1 population (N1: 3.75 and N2: 2.9 days). At the two highest 
copper ion concentrations, the N strains reached their maximum growth earlier than 
the Cu strains (150 uM: N1: 2.41 and N2: 3.31 days; 200 uM: N1: 2.20 and N2: 2.06 
days) (Figure 2 and Table S1). In other words, when exposed to copper ion 
concentration of 150 uM, the N populations exhibited a delay of 0.9 days compared 
with the Cu populations to reach their fastest growth (Table S1). 


3.2. Fy/Fy 


In all cell strains studied, F,\/F,, was negatively affected by increasing copper ion 
concentration (Figure 4), and this effect became more prominent in later days of the 
experiment. As such, the Cu strains showed higher values than the N strains, and the 
strains within each population also differed from one another (as indicated by low 

p values for the effects of the strains in Table 1 and Table 2). The variation in Fy/F), 
was also shaped by two interactions: Experimental Population x Day and Copper 
Concentration x Day (Table 1). When the two population types were analyzed 
separately, the effects were similar, except for the lack of Copper Concentration x 
Day Interaction in the N populations (Table 2). In both cases, the differences 
between the strains remained significant. 
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Figure 2. Growth of Chlamydomonas reinhardtii CW15 cells on media supplemented with different copper ion concentrations 
(0.25-200 uM), measured in terms of the maximal fluorescence of chlorophyll (F,,). N1 and N2 - control strains; Cul and Cu2 - 
copper-adapted strains (four biological replicates). 
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Figure 3 Asymptotic values of the logistic growth curves at different copper ion 
concentrations (0.25-200 uM). N1 and N2 - control strains; Cul and Cu2 - 
copper-adapted strains. 


4. Discussion 
4.1. Adaptation to Copper lons 


Adaptation at the population level can be achieved through altered frequency of 
adaptive alleles, action of gene regulatory networks, or mutations in single elements 
responsible for homeostasis between the cells and environment. There are two 
popular models of adaptation: selective sweep and infinitesimal models. The selective 
sweep model assumes that the changes in allele frequency are rapid and indicative of 
increased fitness of the populations. The infinitesimal model links rapid adaptation 
to many small, parallel changes in allele frequency in response to environmental 
changes (Barghi et al., 2020). 


The toxicity of high copper ion concentrations differs across groups of 
photosynthetic organisms (Beauvais-Fliick et al., 2019). Exposure to high copper ion 
concentration alters ion uptake and intracellular distribution and activates copper 
detoxification mechanisms. For instance, COPT1 and COPT2 are inhibited, while 
other genes responsible for the relocation of copper ion within the cell are activated 
(Andrés-Colas et al., 2006; Sancenén et al., 2004). Recent studies have shown that 
the expression of copper proteins (Cu/Zn superoxide dismutase and plastocyanin) is 
regulated via miRNA synthesis (Burkhead et al., 2009; Pilon, 2017). Similar to other 
noxious factors, copper induces ROS generation, and this response is stronger than 
that of other heavy metals. Contreras-Porcia et al. (2011) noted increased expression 
of antioxidant enzymes (ascorbate peroxidase, glutathione-dependent peroxiredoxin, 
thioredoxin, and glutathione S-transferase) following copper ion exposure. 


Our results, particularly the differences among single cell strains, may be explained 
by the complexity of the copper homeostatic system. Although the specificity of the 
present study did not allow us to highlight specific changes occurring during the 

experiment, we propose possible mechanisms underlying the observed adaptations. 


4.2. Growth Dynamics (F,y) 


Chlorophyll fluorescence measurements have been commonly used since long in 
plant physiology and ecophysiology. They are fast and noninvasive and can provide 
useful information concerning light energy conversion in the photosynthetic 
apparatus (Baker, 2008; Maxwell & Johnson, 2000). 
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Figure 4 Changes in the maximal quantum yield of photosystem II (F\/F,,) from Days 4 to 7 after the inoculation of 
Chlamydomonas reinhardtii CW15 cells on media supplemented with different copper ion concentrations (0.25-200 uM). N1 and 
N2 - control strains; Cul and Cu2 - copper-adapted strains (four biological replicates). 
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Cell strains obtained from the adapted and nonadapted populations differed in 
terms of the maximal fluorescence of chlorophyll a (Figure 2 and Figure 3). 
Differences between the asymptote values of F,,, which could even be observed at 
low copper concentrations, suggest that the selected strains differed in their ability to 
exploit macro- and micronutrients present in the growth medium (Figure 3). 

The more efficient uptake and transformation of nutrients may be a result of changes 
in the synthesis of membrane transporters (Burkhead et al., 2009). We speculate that 
the ability to take up all essential ions from the medium was higher in 
copper-adapted strains than in nonadapted ones. The reverse mechanism, based on 
the reduction of bivalent cations taken up from the environment, limits the ability to 
use other bivalent nutrients (e.g., Ca**, Mg”*, and Zn7*). Other studies have shown 
that at high concentrations (>150 uM), copper ions more strongly inhibited the 
ability of Chlamydomonas reinhardtii cells to take up nitrate from the growth 
medium than Cd** or Zn?* (Devriese et al., 2001). 


The complexity of the copper homeostatic system (uptake, relocation, and 
neutralization) (Burkhead et al., 2009) may be the source of the variability observed 
between the Cu strains. At intermediate copper concentrations, difference in the 
asymptote growth value between strains Cul and Cu2 derived from the same 
population was smaller than that between copper-adapted and nonadapted strains 
(Figure 3). 


Heavy metals can decrease the rate of cell proliferation and induce anomalies during 
cell division, either directly or via ROS generation. According to Luis et al. (2006), 
copper ions interrupt the cell cycle by inhibiting cyclin-dependent kinases. ROS 
interact with nuclear transcription factors to interfere with processes regulated by 
mitogen-activated protein kinases (Leonard et al., 2004). In the present study, 

the cell division rate (estimated based on the initial growth rate) was decreased, 
particularly in strain N1 exposed to copper ion concentrations of 25 and 50 uM. 
Likewise, Jamers et al. (2013) observed a significant decrease in the growth rate 
following exposure to 125 uM CuSO, (but not to <55 uM CuSO,); however, the 
authors did not use a logistic model to better understand the effects of copper ions 
on population growth. In the present study, high copper ion concentrations 
dramatically inhibited the growth of nonadapted strains, which prevented us from 
performing a more detailed analysis of their growth curves. 


4.3. Photosynthetic Performance (F\/F,y) 


Most of the copper ions present in plant cells accumulate in chloroplasts, which is 
directly related to the contribution of these ions to redox processes (Burkhead et al., 
2009). In the view of constant exposure to the harmful effects of copper ions, 
chloroplasts employ various protective mechanisms, including copper binding 

(e.g., metallochaperones) and other protective strategies against the detrimental 
effects of ROS. Specifically, antioxidant enzymes (catalase, peroxidase, and 
dismutase), small-molecule antioxidants (ascorbate, glutathione, preny]l lipids, and 
carotenoids), and enzymes responsible for the regeneration of oxidized antioxidants 
(e.g., monodehydroascorbate reductase and glutathione reductase) enable efficient 
ROS neutralization (e.g., Foyer & Shigeoka, 2011; Nowicka et al., 2016b). 


The thylakoid membrane in the chloroplast, specifically PSII, is the primary site 
affected by copper ion toxicity (Burda et al., 2003, 2006; Yruela, 2009). Therefore, the 
measurement of the maximal quantum yield of PSII based on the fluorescence of 
chlorophyll a (Fy/F\y) is widely used in experiments describing the interactions of 
heavy metals with the photosynthetic apparatus (Luis et al., 2006; Maxwell & 
Johnson, 2000; Nowicka et al., 2016a). Here, the parameters of the fitted growth 
model indicated that copper ion concentration, culture age, and population type, 
as well as the interactions of these factors, affected the photosynthetic apparatus 
(Table 1 and Table 2). The increase in the maximal quantum efficiency of PSII 
suggests that the process of adaptation was continuously underway during algal 
cultivation (factor: day in Table 1). Although mechanisms of chloroplast protection 
(e.g., ROS scavenging) are highly efficient, we cannot certainly pinpoint the precise 
mechanism underlying this adaptation. Notwithstanding, during the experiment, 
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the interactions between copper and PSII evidently differed between the Cu and N 
populations as well as among the strains derived from these populations. 


5. Supplementary Material 


The following supplementary material is available for this article: 


Table S1. Growth parameters [growth rate, inflection point (ie., time point of 
maximum growth), and asymptote] fitted to the strains and presented separately for 
each copper ion concentration. 


Table $2. Additional data to Table 1 in the main text. 


Figure $1. Growth of Chlamydomonas reinhardtii CW15 cells on media 
supplemented with different copper ion concentrations (0.25-200 uM) during the 
first 4 days, measured in terms of the maximal fluorescence of chlorophyll (Fy). 


References 


Andrés-Colas, N., Sancenon, V., Rodriguez-Navarro, S., Mayo, S., Thiele, D. J., Ecker, J. R., 
Puig, S., & Pefiarrubia, L. (2006). The Arabidopsis heavy metal P-type ATPase HMA5 
interacts with metallochaperones and functions in copper detoxification of roots. 

The Plant Journal, 45(2), 225-236. https://doi.org/10.1111/j.1365-313X.2005.02601.x 

Baker, N. R. (2008). Chlorophyll fluorescence: A probe of photosynthesis in vivo. Annual 
Review of Plant Biology, 59, 89-113. 
https://doi.org/10.1146/annurev.arplant.59.032607.092759 

Barghi, N., Hermisson, J., & Schlétterer, C. (2020). Polygenic adaptation: A unifying 
framework to understand positive selection. Nature Reviews Genetics, 21, 769-781. 
https://doi.org/10.1038/s41576-020-0250-z 

Beauvais-Fliick, R., Slaveykova, V. L., & Cosio, C. (2019). Comparative study of Cu uptake and 
early transcriptome responses in the green microalga Chlamydomonas reinhardtii and 
the macrophyte Elodea nuttallii. Environmental Pollution, 250, 331-337. 
https://doi.org/10.1016/j.envpol.2019.04.032 

Burda, K., Kruk, J., Schmid, G. H., & Strzatka, K. (2003). Inhibition of oxygen evolution in 
photosystem II by Cu(II) ions is associated with oxidation of cytochrome b559. 
Biochemical Journal, 371, 597-601. https://doi.org/10.1042/bj20021265 

Burda, K., Kruk, J., Strzatka, K., Stanek, J., Schmid, G. H., & Kruse, O. (2006). Mossbauer 
studies of Cu(II) ions interaction with the non-heme iron and cytochrome b559 in 
Chlamydomonas reinhardtii PSI minus mutant. Acta Physica Polonica A, 109, 237-247. 
https://doi.org/10.1016/S0014-5793(02)03895-4 

Burkhead, J. L., Gogolin Reynolds, K. A., Abdel-Ghany, S. E., Cohu, C. M., & Pilon, M. (2009). 
Copper homeostasis. New Phytologist, 182(4), 799-816. 
https://doi.org/10.1111/j.1469-8137.2009.02846.x 

Contreras-Porcia, L., Dennett, G., Gonzalez, A., Vergara, E., Medina, C., Correa, J. A., & 
Moenne, A. (2011). Identification of copper-induced genes in the marine alga Ulva 
compressa (Chlorophyta). Marine Biotechnology, 13, 544-556. 
https://doi.org/10.1007/s10126-010-9325-8 

Devriese, M., Tsakaloudi, V., Garbayo, I., Leén, R., Vilchez, C., & Vigara, J. (2001). Effect of 
heavy metals on nitrate assimilation in the eukaryotic microalga Chlamydomonas 
reinhardtii. Plant Physiology and Biochemistry, 39(5), 443-448. 
https://doi.org/10.1016/S0981-9428(01)01257-8 

Foyer, C. H., & Shigeoka, S. (2011). Understanding oxidative stress and antioxidant functions 
to enhance photosynthesis. Plant Physiology, 155(1), 93-100. 
https://doi.org/10.1104/pp.110.166181 

Gotelli, N. J. (2008). A primer of ecology. Sinauer Associates. 

Harris, E. H. (2009). The Chlamydomonas sourcebook: Introduction to Chlamydomonas and 
its laboratory use (Vol. 1). Academic Press. 

Jamers, A., Blust, R., De Coen, W., Griffin, J. L., & Jones, O. A. H. (2013). Copper toxicity in 
the microalga Chlamydomonas reinhardtii: An integrated approach. Biometals, 26(5), 
731-740. https://doi.org/10.1007/s10534-013-9648-9 

Jiang, Y., Zhu, Y., Hu, Z., Lei, A., & Wang, J. (2016). Towards elucidation of the toxic 
mechanism of copper on the model green alga Chlamydomonas reinhardtii. 
Ecotoxicology, 25(7), 1417-1425. https://doi.org/10.1007/s10646-016-1692-0 

Leonard, S. S., Harris, G. K., & Shi, X. (2004). Metal-induced oxidative stress and signal 
transduction. Free Radical Biology and Medicine, 37(12), 1921-1942. 
https://doi.org/10.1016/j.freeradbiomed.2004.09.010 


Plucinski et al. / Chlamydomonas — Diversity in Stress Response 


Luis, P., Behnke, K., Toepel, J., & Wilhelm, C. (2006). Parallel analysis of transcript levels and 
physiological key parameters allows the identification of stress phase gene markers in 
Chlamydomonas reinhardtii under copper excess. Plant, Cell & Environment, 29(11), 
2043-2054. https://doi.org/10.1111/j.1365-3040.2006.01579.x 

Maxwell, K., & Johnson, G. N. (2000). Chlorophyll fluorescence - A practical guide. Journal of 
Experimental Botany, 51(345), 659-668. https://doi.org/10.1093/jexbot/51.345.659 

Nowicka, B., Plucinski, B., Kuczynska, P., & Kruk, J. (2016a). Physiological characterization of 
Chlamydomonas reinhardtii acclimated to chronic stress induced by Ag, Cd, Cr, Cu and 
Hg ions. Ecotoxicology and Environmental Safety, 130, 133-145. 
https://doi.org/10.1016/j.ecoenv.2016.04.010 

Nowicka, B., Plucinski, B., Kuczynska, P., & Kruk, J. (2016b). Prenyllipid antioxidants 
participate in response to acute stress induced by heavy metals in green microalga 
Chlamydomonas reinhardtii. Environmental and Experimental Botany, 123, 98-107. 
https://doi.org/10.1016/j.envexpbot.2015.11.008 

Pilon, M. (2017). The copper microRNAs. New Phytologist, 213(3), 1030-1035. 
https://doi.org/10.1111/nph.14244 

Plucinski, B., Waloszek, A., Rutkowska, J., & Strzatka, K. (2018). Copper excess-induced large 
reversible and small irreversible adaptations in a population of Chlamydomonas 
reinhardtii CW15 (Chlorophyta). Acta Societatis Botanicorum Poloniae, 87(1), Article 
3569. https://doi.org/10.5586/asbp.3569 

Sabatini, S. E., Juarez, A. B., Eppis, M. R., Bianchi, L., Luquet, C. M., & Rios de Molina, M. C. 
(2009). Oxidative stress and antioxidant defenses in two green microalgae exposed to 
copper. Ecotoxicology and Environmental Safety, 72(4), 1200-1206. 
https://doi.org/10.1016/j.ecoenv.2009.01.003 

Sancenon, V., Puig, S., Mateu-Andrés, I., Dorcey, E., Thiele, D. J., & Pefiarrubia, L. (2004). 
The Arabidopsis copper transporter COPT1 functions in root elongation and pollen 
development. Journal of Biological Chemistry, 279(15), 15348-15355. 
https://doi.org/10.1074/jbc.M313321200 

Szivak, I., Behra, R., & Sigg, L. (2009). Metal-induced reactive oxygen species production in 
Chlamydomonas reinhardtii (Chlorophyceae). Journal of Phycology, 45(2), 427-435. 
https://doi.org/10.1111/j.1529-8817.2009.00663.x 

Wickham, H. (2016). geplot2: Elegant graphics for data analysis. Springer. 
https://doi.org/10.1007/978-3-319-24277-4 

Yruela, I. (2009). Copper in plants: Acquisition, transport and interactions. Functional Plant 
Biology, 36(5), 409-430. https://doi.org/10.1071/FP08288 


Acta Societatis Botanicorum Poloniae / 2021 / Volume 90 / Article 908 10 


Publisher: Polish Botanical Society 


